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ABSTRACT 

We present photometry of the transient associated with GRB 030329 obtained with 
the CTIO 1.3-meter telescope and the ANDICAM instrument, a dual optical/infrared 
imager with a dichroic centered at one micron. Without the need for light curve in- 
terpolation to produce snapshot broadband spectra, we show that the transient spec- 
trum remained statistically achromatic from day 2.7 to day 5.6, during a re-brightening 
episode. Associating the light in these early epochs with the GRB afterglow, we infer a 
modest level of extinction due to the host galaxy in the line-of-sight toward the GRB: 
Av(h.ost) = 0.30 ± 0.03 mag for (3 = —0.5 and .Ay (host) < 0.4mag (3 a) for any phys- 
ically plausible value of f3 (with flux f v cx A - ^ 3 ). We conclude that the spectral slope 
of the afterglow component was more than (3 = —0.8 between day 2.7-5.6 after the 
GRB, excluding the possibility that the synchrotron cooling break passed through the 
optical/IR bandpass over that period. Taking extinction into account, a decomposition 
of the light curve into an afterglow and supernova component requires the presence 
of a supernova similar to that of SN 1998bw, an afterglow that shows some evidence 
for a second break around day 8-10, and a fifth re-brightening event around day 15. 
Assuming an SN 1988bw-like evolution and a contemporaneous GRB and SN event, the 
peak SN brightness was My = (—19.8 ± 0.4) — 5 log 10 h^ mag. 



Subject headings: gamma rays: bursts — supernovae: individual (1998bw, 2003dh) 
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1. Introduction 

The low redshift long-duration cosmological 7-ray burst, GRB 030329 (Vanderspek et al. 2003) 
afforded an unprecedented glimpse into the aftermath of the explosion. The sheer brightness of the 
early optical afterglow (Peterson & Price 2003; Torii 2003) allowed for high-precision (~ 1 — 5%) 
photometric measurements to be obtained from a bevy of meter-class telescopes (Uemura et al. 
2003; Burenin et al. 2003; Rykoff & Smith 2003; Price et al. 2003a; Lipkin et al. 2003b). The 
optical light curve exhibited 10-30% brightness undulations during the first several hours (Uemura 
et al. 2003) and then showed a prominent break at day 0.48 (Garnavich et al. 2003b; Price et al. 
2003b). Later, the optical transient underwent a brightness resurgence (at about day 2) followed 
by an overall decline that was punctuated by four major re-brightening episodes (Li et al. 2003; 
Lipkin et al. 2003a,c; Sato et al. 2003; Granot et al. 2003). 

Beginning about 6 days after the burst, the first signs of an underlying supernova (SN) compo- 
nent began to emerge in the transient, both spectroscopically (Stanek et al. 2003; Chornock et al. 
2003; Hjorth et al. 2003; Kawabata et al. 2003) and photometrically (Henden et al. 2003; Mathe- 
son et al. 2003); the SN was designated by the IAU as SN 2003dh (Garnavich et al. 2003a). The 
evolution of the SN spectrum in the first two months closely parallels that of the bright type Ic SN 
1998bw (Galama et al. 1998; Patat et al. 2001) thus confirming the previous evidence for (nearly) 
contemporaneous SN and GRB events from the same progenitor in other bursts (see Bloom 2003 
for review). 

While the nature of the SN associated with GRB 030329 is, of course, of intense interest, a 
detailed study of the afterglow holds potential for insight into the structure of GRB explosions. It 
is standard practice to model the afterglow as arising from a jetted point explosion in a constant 
density environment (e.g., Rhoads 1997; Sari et al. 1998) or a wind-stratified circumburst medium 
(e.g., Chevalier & Li 1999). However, the assumption of an instantaneous explosion with a constant 
jet opening angle is untenable both on theoretical grounds and empirical grounds. In particular, 
strong variations have been seen in the early afterglows of GRBs, e.g., GRB 021004 (Fox et al. 2003) 
and GRB 030329 (Price et al. 2003b; Uemura et al. 2003; Burenin et al. 2003; Rykoff & Smith 2003; 
Price et al. 2003a). 

These re-brightenings in GRB 030329 require energy addition to the radiating front. The 
energy addition can be due to slower moving ejecta shells colliding with the earlier emitted but 
faster moving ejecta (Granot et al. 2003); such a hypothesis was first proposed for the variability 
in GRB 021004 (Fox et al. 2003). An alternative interpretation is that the energy addition is 
due to shells with larger opening angles but with slower moving ejecta (Berger et al. 2003). The 
placement of the behavior of the late-time afterglow light curve of GRB 030329 in the context of 
radio modeling can help distinguish between these two hypotheses (Berger et al. 2003). 

An understanding of the multitude of physical processes contributing to the complex behav- 
ior of the transient requires a careful treatment of the high-quality data. Although analysis of 
the vast follow-up data on GRB 030329 in the literature will no doubt be of great use, we have 
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chosen to focus on optical-IR photometric data obtained from the same telescope and reduced 
uniformly. Whereas the burst has been studied extensively at optical wavelengths, only one other 
study presents photometric data in the near- infrared (Matheson et al. 2003). Extending the wave- 
length coverage by a factor of two allows stronger constraints to be placed on the extinction, and 
a better determination of the slope of the afterglow spectrum as a function of time. 

In §2, we present the optical and infrared data obtained with ANDICAM between day 2 and 
day 23 after the GRB trigger. To disentangle the various physical components of the transient 
light, in §3 we first show that while the overall source brightness fluctuated in the first 6 days, the 
spectrum of the transient remained fixed. In §3.1, we then fit this early spectrum to find a constraint 
on the line-of-sight extinction to the GRB and the intrinsic afterglow spectrum. Using these values, 
we then decompose the afterglow and the supernova light in §3.2. We compare our results to those 
obtained in other studies and conclude with some implications for the GRB-supernova connection 
and the explosion structure of GRB 030329. 

2. Observations, Reductions, and Calibrations 

All the observations reported herein were carried out with the ANDICAM instrument mounted 
on the 1.3-m telescope at Cerro Tololo Interamerican Observatory (CTIO). The ANDICAM is a 
dual-channel camera constructed by the Ohio State University instrument group 1 . The camera con- 
tains a dichroic that enables two imagers to be simultaneously illuminated: a Rockwell 1024x1024 
HgCdTe "Hawaii" array, and a Fairchild 447 2048x2048 optical CCD. This instrument was for- 
merly mounted on the Yale 1-m telescope at CTIO, operated by the YALO consortium (Bailyn 
et al. 1999). ANDICAM has recently been transferred to the 1.3-m telescope at CTIO, formerly 
used for the 2MASS survey, where it has been operating since February 2003 by the SMARTS 
consortium 2 . The image scale is smaller at the 1.3-m than it was at the 1-m, and the field of view is 
now 2.4' x 2.4' for the IR array, and 6.3' x 6.3' for the CCD. The optical images are routinely double 
binned to provide a pixel scale of 0.37"/pixel, while the IR channel significantly oversamples the 
seeing (typically 0.8") with 0.14" pixels. 

The field seen by the IR array can be repositioned slightly by adjusting three tilt axes of 
an internal mirror. This allows us to "dither" the IR position while an optical integration is 
underway. These dithered IR images can be used to generate a sky image without actually moving 
the telescope. Such sky images are slightly inferior to those made by moving the telescope, since 
the field is viewed through a different part of the filter, leading to decreased accuracy in flat 
fielding. This problem is most significant in the X-band, which, at the 1.3-m, also suffers from high 
background levels. 



1 See http : //www. astronomy . ohio-state . edu/ANDICAM. 
2 See http : //www. astro .yale . edu/ smarts. 
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Table 1 gives a log of the observations with ANDICAM. The data have been grouped into 
13 individual epochs. As can be seen the total duration of the epochs was always less than 40 
minutes, even over epochs where BVIJH images were acquired. Individual optical images were 
reduced using IRAF/CCDRED 3 . The bias was removed using a mean of several bias frames 
and the images then flat-fielded using the best composite twilight flats for the appropriate filter. 
Twilight flats were not obtained on every night of the transient observations but the differential 
gain variations over several days should not affect the photometric uncertainty at significantly more 
than the 1% level. Images where the OT light was contaminated with cosmic rays were excluded 
from the final sample. 

Individual infrared images were reduced in the following way. Flat fields were created in J- 
and //-bands by subtracting pairs of twilight flats with different sky brightnesses. After dark 
subtraction each science frame was divided by the (normalized) flat field. Next, a sky frame was 
created by taking the median of all normalized images for a given epoch, and subtracted after 
rescaling to match the counts in the original images. In the H- and i^-bands strong "ripples" are 
present in the reduced images. These features likely result from the ANDICAM dithering mode: 
as mentioned, at each mirror position the light passes through a different part of the filter, creating 
residual flat fielding errors. In the H-hanA the features could be successfully removed by creating 
residual flat fields, created from all images (at all epochs) taken at the same mirror position. In 
the if-band the residual pattern dominated the flux of the transient even after this correction; 
therefore, the K-band observations are not considered further in this paper. Finally, for both the 
optical and IR images, cosmic-rays were removed using a Laplacian edge-detection algorithm (van 
Dokkum 2001). 

Usually only one or two images per epoch were obtained per optical filter, while each epoch 
generated between 7-21 dithered images per IR filter. In a given IR image, often only a few objects 
were detected at better than 3 a, requiring a complicated image alignment procedure in order to 
stack the images. In a custom Pyraf code, we first removed any large scale gradient from a copy 
of all the images using IRAF/MKSKYFIT, masked hot and cold pixels, and then smoothed these 
images with a Gaussian of width 3 pixels. The pixel-to-pixel rms was then determined using an 
iterative sigma-clipping algorithm. This rms was then used as an input to DITHER/PRECOR 
to mask out all regions of the sky where less than 18 out of the neighboring 5x5 pixels were above 
three times the rms. The effect was to create masked images where only objects appear, and the 
remainder of the pixels are set to zero. 

An initial estimate of the true offset was computed using the positional and mirror keywords 
in the FITS headers. For epochs where the seven-point dither pattern was repeated at least once, 



3 See http://iraf.noao.edu. IRAF is distributed by the National Optical Astronomy Observatories, which are 
operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

4 See http://www. stsci . edu/resources/sof tware_hardware/pyraf . 
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we found that a given set of three mirror tilts (recorded in the FITS headers) consistently mapped 
source positions to better than one pixel rms. From a training set, we fit a linear least-squares 
transfer (3x2) matrix to map the three tilts to a nominal x,y source position taking into account 
the reported telescope pointing, plate orientation, and pixel scale; we assumed no differential sky 
rotation from image to image. The difference between the nominal x, y position in the fiducial 
image and the nominal x,y position in a given image yielded an estimate of the nominal offset. 

To find the precise offsets between a fiducial image and the other images in a given epoch, we 
cross-correlated the masked images using DITHER/CROSSDRIZZLE. The highest peak in the 
cross correlation within a 50 pixel width box about the nominal offset position was then found and 
centroided using SHIFTFIND; this process yielded a list of sub-pixel offsets with uncertainties less 
than 0.1 pixels. Individual epochs were then stacked with DITHER/DRIZZLE using these offsets 
without re-sampling the data to a finer grid (i.e., drizzle . scale = 1 and drizzle .pixfrac = 
1). The offsets between the stacked images of all IR epochs were then computed using the same 
cross-correlation technique. The final epoch stacks were made using these epoch-to-epoch offsets 
as input to the secondary geometry parameters for drizzle (drizzle . dr2gpar), repeating the entire 
process described above. The result was a set of stacked, registered images for all IR epochs. 

Using IRAF/FITPARAM, we calibrated our first epoch (t = 2.67 day) of optical images to 
secondary field standards reported in Henden (2003). An average extinction vs. airmass curve from 
nearby Paranal 5 was assumed and the differential color terms found to be consistent with zero. 
This implies that our cumulative filter+telescope response curves reasonably approximate those 
of the Henden system, and thus closely resemble the Johnson BV and Kron-Cousins / filters (i.e., 
Landolt system). The photometry for all subsequent epochs were computed differentially from the 
zeropoints established on the first epoch. As a check, primary optical Landolt standard stars were 
observed on several different photometric nights and found to yield consistent photometry to better 
than 5% accuracy. 

The infrared standard star SJ9144 (Persson et al. 1998) was observed on the third (photomet- 
ric) epoch (April 3 UT). Using curve-of-growth aperture photometry, we determined a zeropoint 
for the IR science observations at that epoch and assume a systematic zeropoint uncertainty of 
0.05 mag. Aperture photometry corrected to the "infinite aperture" magnitudes were computed for 
the transient and all detectable secondary standards in the field for all epochs. As with the optical 
observations, photometry was computed differentially from this fiducial epoch. The resulting pho- 
tometry is given in Table 1. Note that while our J-band photometry agrees with the results from 
Matheson et al. to within the statistical uncertainties, as brought to our attention by K. Krisciunas, 
our H-b&nd measurements appear to be about 0.2 mag fainter at similar epochs to the results in 
Matheson et al. Unfortunately, there are no intercomparison stars at common epochs between the 
two datasets. Therefore we were unable to determine the origin of the difference in the .ff-band. 



5 http: //www.eso . org/instruments/isaac/imaging_standards .html. 
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Transient fluxes were first computed from the observed magnitudes assuming the transfor- 
mations given in Fukugita et al. (1995) (optical) and Bessell & Brett (1988) (IR). This revealed 
a small and consistent curvature in the broadband spectrum in the first several epochs, with an 
incident spectrum ranging from (3 ~ —1 at the blue end to (3 ~ —0.6 in the red end. Since the 
effective wavelength and zeropoint of filters depend upon the incident spectrum we recomputed 
these quantities by fitting a 2nd-order polynomial to the first epoch. This fiducial spectrum was 
then input to IRAF/CALCPHOT, part of the STSDAS/SYNPHOT package 6 to compute the 
effective wavelengths (A e gf) and Vega-referenced zeropoints of the filters. The difference between 
fiducial filter wavelengths and the effective filter wavelengths was small, a few tens of Angstroms. 
Though the input spectrum was seen to evolve somewhat over the course of our observations (see 
below), given the small changes in A e fr from the fiducial values, in the subsequent analysis we fix 
A e g for all filters at all epochs. The resulting transient fluxes are given in the last column of Table 
1. 



3. Results 

With dense high-precision sampling, the complexity of the early light curve of GRB 030329 
was revealed. The temporal features of the early afterglow (t < 1 day) have been described in 
detail elsewhere (e.g., Uemura et al. 2003; Lipkin et al. 2003a). Later time "bumps and wiggles" 
(or rebrightenings) are also detected in the ANDICAM light curves, shown in Figure 1. 

Figure 2 shows the transient broadband spectra over the 13 epochs of observations, uncorrected 
for the effects of dust reddening. As evident, while the overall brightness fluctuated (Fig. 1), there 
was no apparent change in the spectrum over the first several epochs. To understand the statistical 
significance of this apparent achromaticity, we take an empirical approach by examining the changes 
in the 10 color indices, constructed from the 5 observing filters (see Fig. 6). We compute the color 
change from the first epoch for all subsequent epochs and estimate the significance of change 
from the statistical errors on the measurements (i.e., only the differential errors are considered). 
Between epochs 1 through 4, 28 out of 30 color changes are consistent within ±2 a from zero 
change; the significance of the B — V and B — I colors changes appears high on epoch 3 (by 3 and 
4d, respectively). From epoch 5 onward, the distribution of color change significances fans out (as 
will be shown in §3.2, there are also secular trends in the evolution of some colors). Specifically, 
the Kolmogorov-Smirnov (KS) probability that the distribution of color change significances from 
epochs 1-4 is drawn from a Gaussian of sigma unity centered about zero is P«s = 0.5 but then 
plummets to below 0.05 by day 10. We therefore conclude that the transient spectrum did not 
begin to evolve significantly until after the fourth epoch of observations. 

The lack of a significant color change in the early epochs allows us to decouple the analysis of 



See http : //stsdas . stsci . edu/STSDAS . html. 
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the afterglow (dominating the first four epochs) and the supernova (contaminating the remainder of 
the epochs). In what follows, using the first four epochs only, we determine the intrinsic afterglow 
spectrum and the line-of-sight extinction concurrently. With these results, we then decompose the 
afterglow component and the supernova component. In accordance with early spectroscopic studies 
(e.g., Matheson et al. 2003), we find that any contribution to the total optical/IR light from an 
underlying supernova was < 5% in all filters during the first four epochs (see §3.2). 

3.1. Constraints on the Line of Sight Extinction 

The achromaticity of the transient spectrum is consistent with the expectation of afterglow 
emission in the standard synchrotron shock model (e.g., Sari et al. 1998). In the model, the intrinsic 
optical-IR afterglow spectrum is dominated 7 by a featureless power-law spectrum characterized as 
f u oc A _/3 . There is, however, an apparent downward curvature to the observed broadband afterglow 
spectrum (Fig. 2), which can modeled in the first four epochs as f v « (—3.17 ± 0.21)A3 + (119 ± 
4)A 3 - (131 ± 17) juJy, from A » 4000 A to 18000 A (A 3 = A/1000 A). 

Given the expectation of an underlying power-law, the curvature is likely due to extinction by 
dust along the line-of-sight to the burst. We constrain the intrinsic extinction by normalizing the 
first four epochs to a unity flux density in /-band and fit the function, 

f v (A) (observed) = C\~ p x T [Ry (Gal) , Ay (Gal) , A] x T[i?y(host), Ay (host), A/(l + z)\ /xJy, (1) 

where Ry and Ay are the extinction curves and equivalent V-band extinction through the Galaxy 
("Gal") and the line-of-sight through the host galaxy/progenitor system ("host"), respectively. 
The redshift of the host is z = 0.1685 (Greiner et al. 2003). A parametrized transmission function 
T[Ry, Ay, A] is adopted from Cardelli et al. (1989). 

The Galactic extinction toward the transient, estimated from Galactic foreground IR emission, 
is E(B — V) = 0.025 mag (Schlegel et al. 1998). Following from a recent analysis by Burstein (2003), 
the statistical error on this estimate is (Te(b-v) ~ 0.002 mag, or about 10%. In the extinction fit 
we fix fly (Gal) = 3.1 and thus fix Ay (Gal) = 0.083 ± 0.017 mag. 

By minimizing x 2 , we first fit simultaneously for the values of C, Ry (host), Ay (host), and 
(3, and found the preferred value of Ry (host) is small (< 2); we consider this unlikely to be the 
true value of Ry given the indication (from spectroscopy and imaging) that the host is a starburst 
dwarf galaxy (e.g., Matheson et al. 2003): for such galaxies we would expect an extinction curve 
Ry ^ 3. Therefore, we examined the results for two values of -Ry(host) = 3.1 and 5.1 (LMC-like). 



7 Though the early spectra show atmospheric absorption and star-forming emission lines from the host (Stanek 
et al. 2003; Hjorth et al. 2003), these narrow features are unlikely to dominate the wide filter light, as the equivalent 
width of the lines (<; 50 A) is significantly smaller than the width of the filters (j> fOOO A) 
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For i?y (host) = 3.1 (5.1), the best fit values are Ay (host) = 0.94 ± 0.24 mag (1.8 ± 1.0 mag), 
(3 = 0.11 ± 0.22 (0.5 ± 0.8), * 2 /dof = 0.44 (1.34). 

Note that the errors given are only derived from the diagonal elements of the cross-correlation 
matrix. As Figure 3 demonstrates, there is a strong covariance between Ay (host) and (3 in the 
sense that an intrinsically more blue afterglow requires a higher value of extinction. 

Though the fits given above are formally acceptable, in both cases, the intrinsic afterglow 
spectrum would be unusually blue if the origin is a relativistic synchrotron shock. In previously 
studied afterglows, the observed temporal decline index a (f u oc t a ) could be used to constrain the 
intrinsic spectral index through the a— (3 closure relations prescribed from synchrotron shock physics 
(see Price et al. 2002 and references therein). In GRB 030329, however, the complex temporal 
behavior does not yield an obvious estimation for the value of [3. That said, if we adopt the 
interpretation of the light curves from Granot et al. (2003), then the temporal break observed at 
day 0.48 ± 0.03 (Price et al. 2003b) can be attributed to a jet break in the standard manner, and 
(3 then equals — (a.2 — l)/2. Since Price et al. (2003b) find the decay after the break to be oc2 ~ 2, 
we expect that (3 « —0.5. 

A similar value (an upper limit actually) for (3 is inferred by a consideration of the range of 
values of the electron energy index spectrum p, where the number density of shocked electrons as a 
function of energy E is n(E) cx E~ p . For the energy in the electrons to be finite 8 , p must be greater 
than 2. In the simplest homogeneous density model where v c > u(opt), (3 < —0.5, consistent with 
value inferred above from the temporal break. The cooling frequency v c corresponds to an energy 
at which electrons have radiated a significant fraction of their initial energy in the lifetime of the 
shock (e.g., Sari et al. 1998). We also note that the apparent achromaticity of the afterglow over 
the first 4 epochs, implies that no synchrotron break frequency (see, e.g., Sari et al. 1998) moved 
through the optical/IR bandpass from 2.7 to 5.6 days after the burst. The cooling break frequency, 
then, likely remained above ~ 10 14 Hz during this time. Imposing (3 < —0.5 and -Ry(host) = 3.1, 
we find that Ay(host) < (0.30 ± 0.03) mag. If we consider any fit acceptable with x 2 /dof < 2, 
following from Figure 3, the extinction from the host can be zero for p ~ 2.6. A similar value for 
the extinction constraint — Ay (host) < (0.36 ±0.04) mag — is found for Ry (host) = 5.1. If (3 = —0.5 
(see Fig. 4), the combined transmission of the Galactic and extragalactic dust is: Tb = 0.60, 
T v = 0.68, Tj = 0.77, Tj = 0.89, and T H = 0.92. 

Tiengo et al. (2003) and Lipkin et al. (2003c) have noted that the X-ray to optical ratio of the 
afterglow remained approximately constant from 5 hours to 60 days after the burst, with a spectral 
index (if indeed connected by a single power-law) of O pt-x ~ — 1. This single power-law is clearly 
inconsistent with our assertion that the optical-IR spectral slope is greater than —0.8. This suggests 
that a (cooling) break is required between the optical and X-ray bands above which (3 < — 1 (see, 



8 A value of p < 2 would not violate the finite energy requirement as long as there exists a high-energy cut-off in 
the electron energy spectrum (Dai & Cheng 2001). 
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for example, Fig. 2 of Bloom et al. 1998). This statement is at odds with the conclusions of Tiengo 
et al. (2003). If indeed, the X-ray to optical flux ratio remained fixed during our first four epochs 
(Lipkin et al. 2003b), this implies that v c remained constant over that period. In turn, this lends 
support to the hypothesis (Berger et al. 2003) that a jet break — from the fast-moving ejecta 
responsible for optical/IR and X-ray emission — occurred before this time. This statement is at 
odds with the conclusion of Lipkin et al. (2003c). 

3.2. Decomposing the Afterglow and the Supernova 

With an empirical examination of the color evolution in the transient, we showed that the 
spectrum began to change after day ~6. We can now confirm this with a physical model for the 
transient, assuming a single power-law afterglow spectrum and dust extinction. Fixing the values 
of Tj found above, we fit (3 as a function of time. The result is shown in Figure 5. As seen, the 
first 5 epochs are consistent with a single power-law (3 = —0.5, synchrotron emission from the 
GRB afterglow. Then, the transient spectrum becomes inconsistent with a single power-law as, 
presumably, the light from the underlying supernova begins to contaminate. 

While the spectrum on epoch 5 is formally consistent (at the 1.5 a level) with (3 = —0.5, 
Figure 2 shows that there was a small apparent brightening in the V- and .ff -bands relative to 
the spectrum index between the I- and J-bands. The -ff-band brightening could be the same 
"color event" described in Matheson et al. (2003); the brighter l/-band flux could be accounted for 
by the onset of a non-negligible (~5— 10%) contribution from the associated SN. In the subsequent 
epochs, the broadband spectrum clearly evolves, with the light from the SN becoming an increasing 
fraction of the total light at optical wavelengths (see Matheson et al. 2003 for details of the relative 
contribution between the SN and the afterglow derived from spectroscopy). 

To decompose the afterglow and supernova light, we constructed a set of predicted SN light 
curves in the BVI bandpasses by dimming the magnitudes of SN 1998bw to the redshift of 
GRB 030329 (Galama et al. 1998; Patat et al. 2001). We calculated k-corrections derived from 
1998bw spectra by using the prescription as described in Kim et al. (1996), calculating the k- 
correction at each 1998bw spectral observation epoch. A cumulative line-of-sight extinction to 
1998bw of E(B — V) = 0.061 was assumed (Woosley et al. 1999). To produce a set of synthetic SNe 
light curves at the (time-dilated) epoch of our GRB 030329 observations, these individual 1998bw 
epochs were fit with a least-squares spline, with the uncertainty at any epoch gauged from the rms 
scatter about the fit. The calculated k-corrections transform the restframe (1998bw) bandpasses 
of U, B, V, R, I to the observed bandpasses (GRB 030329) of B, V, R, I, Z, respectively. The 
predicted magnitudes versus time were converted to flux [/sn(Aj, ij); units of /uJy]. To predict the 
IR magnitudes, we fit a time-evolving blackbody spectrum to the synthetic SN optical light curves; 
this temperature showed a reasonable hard-to-soft evolution from log T(K°) = 6.1 to logT(K°) = 
5.8. We then converted the blackbody fluxes to magnitudes, giving a crude estimate of the J- and 
//-band magnitudes of the SN component. 
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For a fixed dimensionless normalization of the SN model, csn, we fit for the scaling of the 
afterglow component, c Jia f t , at each epoch j where the observed flux in filter i is, 

, 0.5" 

CSN /sN(Aj, tj) + Cj ja ft(tj) 



f v (\ i ,t j )\jiJy] = T i x 



eff,I— band 



(2) 



The first term is a scaled model of the supernova component, with csn constant for all epochs. The 
second term is the afterglow component, referenced to the flux at the effective wavelength of the 
/-band filter (X e s,i -band'^ see Table 1). The value Cj^(tj) (units of /txJy) varies as a function of 
time, but the afterglow spectrum is fixed with (3 = —0.5. The filter transmission values are given in 
§3.1. Since pre-imaging of the optical field of GRB 030329 suggested a faint host galaxy compared 
with the observed transient magnitudes (i? n0 st > 22.5 mag; Wood-Vasey et al. 2003; Blake & Bloom 
2003) we do not include a constant host component in the fit. 

Figure 6 shows the observed and modeled color evolution. Without an SN component, the 
"afterglow only" model cannot account for the secular evolution in colors, which become manifest 
after about day 6. Instead, the addition of a 1998bw-like SN — brighter by a factor of csn ~ 1.5 
— appears to better represent the observed color evolution, notably from day 6 to day 13. A value 
of csn = 3 clearly over-predicts the color changes. Had a broadband color change been due to a 
change in the afterglow synchrotron spectrum itself, all color indicies should have shown a near 
simultaneous change; for the passage of a cooling break, all colors would become more red. Instead, 
some colors become more red while others become more blue, consistent with a SN source spectrum 
which peaks in the l/-band. 

In Figure 7, we show the fit values of Cj^a for three values of csn as well as the SN model. 
This analysis demonstrates that the contribution of the SN was <; 5% at all the epochs used to 
derive the extinction in §3.1, and hence the determination of the values for Ay and (5 are not 
strongly affected by the SN. In all cases, the fourth re-brightening event of the afterglow (at t ~ 5.2 
day; bump "D" in Granot et al. 2003) is clearly seen. With a non-negligible contribution from the 
SN, the decomposition reveals another steepening around 8-10 days, followed by a possible fifth 
re-brightening event around 15-20 days. Our favored valued of csn = 1.5 implies an equipartition 
of the afterglow and supernova contribution to the transient light at day 11 — 13 in the /-band, 
consistent with the spectral decomposition from Matheson et al. (2003). 



4. Discussion and Conclusions 

Our analysis of the ANDICAM data has shown that, until day ~5, the transient light is dom- 
inated by the afterglow. During the early epochs, the spectrum remains statistically achromatic, 
even during a re-brightening event. Starting around day 6, the transient slowly evolves in color as 
the supernova begins to contribute to the total light (Fig. 6 and §3.2). In the optical bandpasses, 
the early broadband spectrum is consistent with (3 = —0.94, reported by Stanek et al. (2003) and 
Lee et al. (2003) but is more shallow in the IR (Fig. 2). Hjorth et al. (2003) report a significantly 
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more red spectral index of the afterglow component ((3 = —1.2 ± 0.05) constrained mostly from 
their April 3.10 spectrum obtained on the Very Large Telescope (VLT). This is inconsistent with 
our broadband spectral measurement. One explanation for the discrepancy is that there may have 
been considerable differential light loss (Filippenko 1982) in the VLT spectrum: over the course of 
the three 600 sec exposures 9 , the parallactic angle (East of North) changed from ss —170° to 170° 
yet the position angle of the slit was fixed at 123.6°, i.e,. more than 40° from the optimal angle. 
Though the spectroscopic slit width was between 1.3 to 2 times the average seeing, the combination 
of the high airmass (secz = 1.44) and this chosen slit angle could easily account for the anomalous 
redness of the Hjorth et al. (2003) spectra. 

The measurement of the line-of-sight dust extinction in early GRB afterglows is critical to 
determine the parameters of the shock physics (e.g., Berger et al. 2001), to provide a means to 
deredden the flux of any coincident supernova (e.g., Price et al. 2002), and (in conjunction with 
an Nh measurement from X-ray spectroscopy) to constrain the extent of dust destruction near 
the explosion site (Galama & Wijers 2001). Our constraint on the line-of-sight extinction toward 
GRB 030329, Ay(host) < (0.30 ± 0.03) mag is obtained under the assumption that the electron 
energy index (3 < —0.5. Matheson et al. (2003), using independent broadband photometry from 
U-H-b&nd on day 5.6, find a consistent value of Ay (host) = 0.23 ± 0.25 and (3 = —0.71 ± 0.21. As 
demonstrated in Figure 3, the coupling between the fit values of Ay (host) and (3 is non- negligible; 
therefore, we caution that any subsequent use of these measured quantities should take this covari- 
ance into account. Specifically, the lowest values of the x 2 surface are well-reproduced under the 
simple constraint that (3 = 0.92 Ay (host)-0. 75 (for _Ry(host) = 3.1). Thus if Ay (host)= Omag, 
then (3 = -0.75. 

We compare the line-of-sight extinction toward GRB 030329 with the average extinction to- 
ward Hn regions in the host galaxy. Assuming case-B recombination and the Osterbrock (1974) 
interstellar extinction curve, the Ha/H/3 ratio can be converted to E(B — V) using 



Using Balmer line luminosities measured by (Hjorth et al. 2003) we find E(B — V) = —0.01 ± 0.25. 
Taking the Galactic extinction of E(B — V) = 0.025 into account we find E(B — V) < 0.23, and 
Ay < 0.76. This is comparable to the extinction inferred by Matheson et al. (2003). We note that 
the line luminosities have not been corrected for Balmer absorption, as the host luminosity is not 
yet known; however, since the equivalent widths of the Balmer lines are large when the supernova 
and afterglow still dominated, this correction will be negligible for Ha and H/3. The low extinction 
in the line-of-sight toward GRB 030329 is therefore typical for star forming regions in this galaxy, 
and should therefore not be taken as evidence for the destruction of dust by the GRB itself (e.g., 
Waxman & Draine 2000). 



9 See http : // archive . eso . org. 




(3) 
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As shown in Figure 3, a value of (3 < —0.8 is formally excluded from the locus of acceptable fits 
in both cases of Ry (host). This then excludes any possibility that the synchrotron cooling break 
moved through the optical/IR bandpass from 2-6 days after the GRB. In the §3.2 we showed that 
the passage of cooling break could not explain the subsequent color changes from day 6-15, and 
was therefore also unlikely to have moved through the optical/IR bandpass during that time. 

The suggested late-time afterglow features seen in Figure 7 are also an important clue for 
understanding the structure of GRB explosions. Berger et al. (2003) suggest that the break in the 
radio light curve around day 10 should have been accompanied by a second break in the optical 
light curve. Our decomposed afterglow light curve does indeed show some evidence for such a break 
around day 8-10. The precise timing of the break and the slope of the post-break afterglow light 
curve depends sensitively on the value of csn- Based on only a few later-time data points, we also 
suggest tentative evidence for a fifth re-brightening event around day 15-20. If true, the mechanism 
for the later-time injection of energy into the radiating front must continue at least until this time. 

As shown in Figure 6, a concurrent supernova about 50% brighter than SN 1998bw helps to 
reproduce the reddening in the B — V color and the bluing in the V — I and V — J colors. The 
supernova plus afterglow model is certainly not a perfect (nor unique) fit to the data, likely due 
to the inherent uncertainties in modeling the supernova a priori. For instance, the modeled SN 
over-predicts the late-time flux at epoch 13 by ~ 0.2 — 0.3 mag (this effect was also noted by Hjorth 
et al. 2003). While the spectroscopy of the SN shows it to be conclusively of type Ic (Hjorth 
et al. 2003; Chornock et al. 2003; Matheson et al. 2003; Kawabata et al. 2003), there is a large 
observed diversity in color and light curves from this class of supernovae (e.g., Iwamoto et al. 1998; 
Mazzali et al. 2002). SN 2003dh could simply have risen and decayed faster than 1998bw (see also 
Bloom et al. 2002). Instead, the supernova could have occurred a few days before the GRB (the 
"supranova" scenario; Vietri & Stella 1998; Hjorth et al. 2003). However, as Guetta & Granot 
(2003) emphasize, had a GRB occurred within a few days (or even a few months) after an SN, the 
optical depth to Thompson scattering from the SN shell would be much larger than unity, thus 
extinguishing the GRB. A highly asymmetric SN explosion could produce lower optical depths 
for certain viewing geometries, but the lack of strong polarization in the SN spectrum suggests a 
modest level of asymmetry (Kawabata et al. 2003). We therefore conclude that the SN and GRB 
explosions were contemporaneous to within the free-fall collapse time of the progenitor star (see 
Guetta & Granot 2003). 

Since we did not observe during the peak of the SN (see Figure 1), we do not have a direct 
measurement of the peak SN brightness. However, assuming a 1998bw-like evolution of the super- 
nova component, scaled to be 1.5 times brighter (0.44 mag), the peak absolute magnitude of 2003dh 
is My = —19.8 — 5 log 10 /i65 mag; we estimate an uncertainty in this value due to csn of 0.3 mag. 
There is an additional ~10% systematic uncertainty in this quantity owing to the unknown peculiar 
velocity of the host of SN 1998bw. 

We end by highlighting one of the unique features of ANDICAM, namely the ability to si- 



-13- 



multaneously image afterglows in both the optical and infrared bands. Aside from allowing for 
efficient characterization of the broadband spectrum to measure line-of-sight extinction, any short 
timescale (^minutes) color changes in the afterglow could be detected unambiguously. If such 
observations were to be conducted in the first few hours after a GRB, strong constraints on the 
passage of the synchrotron peak frequency could be obtained, leading to precision measurements 
of the peak Lorentz factor of the shell. In addition, the constraints on short term color variations 
would lead to constraints on the patchiness of clouds near the progenitor (Kumar & Piran 2000). 
Indeed, with such instruments in the Swift era we look forward to a level of insight even beyond 
that gleaned from observations of the remarkable transient of GRB 030329. 
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Time since GRB 030329 [days] 

Fig. 1. — The observed light curve of the transient associated with GRB 030329. The complex 
bumps and wiggles in the early afterglow, especially the re-brightening event at day ~5.2, are 
apparent. Shown are the modeled light curves in the B VI filters (dashed, dot-dashed, dash-dotted, 
respectively) of SN 1998bw at the redshift of GRB 030329 and dimmed by the total line-of-sight 
extinction inferred for the afterglow (see text). The SN brightness scaling (csn) is later varied in 
§3.2 to decompose the light curve into an afterglow and SN component. 
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Fig. 2. — The broadband spectral evolution of the transient associated with GRB 030329. The time since the 
GRB is noted beside each spectrum. The data have not been corrected for reddening. No significant spectral 
evolution was detected in our first four epochs, at least until 6 days after the GRB. For reference, the initial 
optical spectral index reported by others (/3 = —0.94) is shown and is consistent with our observed spectra 
in the bluer filters. However, the IR data show the spectrum to flatten toward longer wavelengths, indicative 
of extinction due to dust reddening. A systematic uncertainty of 4% in the conversion of magnitude to flux 
has been included. 
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Fig. 3. — Coupling between the host line-of-sight restframe extinction measurement and the intrinsic spectral 
power-law index /3 for two different extinction curves [top: i?v(host) = 3.1 (Galactic); bottom: i?v(host)= 5.1 
(LMC)]. The first four epochs are used in this fit, scaled to the same flux value at the V-band. A Galactic extinction 
of j4v(Gal)= 0.083 mag, i?v(Gal) = 3.1 is assumed. Lower values of host extinction spectra result in lower (more red) 
values of (3. The data prefer a blue afterglow with high host extinction, though physical constraints require a more 
red afterglow: the horizontal dashed lines show the predicted /3 for p = 2 in the fast and slow cooling regimes. For 
p > 2, in either regime, (5 must be more negative. Solid contours show the 1-, 2-, 3-, and 5-<r reflect the 2-parameter 
demarcations from the \ surface. The dash-dotted line shows an acceptable region of the data, where xVdof < 2. 
Therefore, /3 must be greater than —0.8 in the first four epochs, ruling out the hypothesis that the cooling break, u c , 
passed through the optical/IR bandpass between days 2 and 6. 
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Fig. 4. — Best-fit afterglow spectrum plus extinction assuming j3 = —0.5 (p = 2, slow cooling). As in 
Figure 3, the Galactic extinction is fixed. The dotted line shows the unreddened afterglow spectrum. 
The dashed line shows the afterglow spectrum reddened only with Galactic extinction. The solid 
line shows the afterglow after reddening from the Galaxy and a host line-of-sight extinction of 
Ay = 0.30 ± 0.03 mag (assuming itV(host) = 3.1). This fit is acceptable, with x 2 /dof = 0.94. 
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Fig. 5. — Power-law fit to the evolving intrinsic transient spectrum as a function of time, fixing 
the extinction measured in Fig. 4. The top panel gives the value for (5 and the associated 1 o 
error assuming an acceptable fit. The bottom panel gives the x 2 /dof of the fit. Measurements 
at epochs marked with solid circles are acceptable fits (x 2 /dof< 2), implying that the transient 
can be adequately modeled as a single power-law. Open circles show unacceptable fits to a single 
power-law. 
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Fig. 6. — Color evolution of the observed transient compared with a model for the afterglow and the supernova 
component. The observed data, uncorrected for dust extinction, is shown with l-cr error bars. The model, using a 
mix of an afterglow and a 1998bw-like supernova, was constructed following §3.2. The dashed horizontal line derives 
from the best fit result to Cj, a ft assuming no contribution from a supernova (i.e., csn = 0). The solid (dotted) lines 
show the color results for csn = 1-5 (3.0). The csn = 1.5 helps reproduce the apparent bluing trend in V — I, V — J, 
and V — H and the reddening trend in B — V . Such color trends are a natural consequence of an SN source spectrum 
which peaks around 5000-6000 A. 
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Fig. 7. — Evolution of the afterglow and supernova brightness at the effective wavelength of the 
/-band filter for different levels of contribution from an underlying supernova. The temporal de- 
pendence on the scaling parameter Cj A ft is shown as a dotted curve assuming no contribution from 
a supernovae (csn = 0); in Figure 6 we show that a model with csn = does not reproduce the 
color evolution adequately. Instead, a reasonable agreement in both color and flux is obtained by 
setting csn = 1.5. The sum of both components for a given csn should equal that of the dotted 
curve within the errors. The fourth re-brightening event ("D" from Granot et al. 2003) is labeled 
with an arrow at day 5.2. There is some evidence for a break in the light curve between day ~8-10 
and another re-brightening event in the light curve around day 15-20, labeled as "E". The fluxes 
have been corrected for the effects of dust extinction. 
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Table 1. Photometric Observations of GRB 030329 



Date a 


Ai b 


Filter 


Exp. Time 


Airmass 


Mag. c 


Flux d 


UT 


day 




sec 


sec z 




fijy 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 










Epoch 1, St — 


38.5 min 






Apr 1 03:33:29 


2.664 


B 


600.00 


1.64 


17.67 ±0.01 


341.1 ±9.9 


Apr 1 03:45:52 


2.673 


V 


600.00 


1.67 


17.26 ±0.01 


442.8 ± 17.7 


Apr 1 04:01:30 


2.684 


I 


600.00 


1.71 


16.48 ±0.01 


613.8 ±31.9 


Apr 1 yJ6.66.Zi 


Z.004 


T 

J 


cfin i n 


l.u4 


10.04 zh U.lz 


ODO.O ± 1U4.1 


Apr 1 03:45:50 


2.673 


H 


560.18 


1.67 


15.02 ±0.05 


981.2 ±62.0 








Epoch 2, 8t = 


35.1 min 






Apr 2 02:19:40 


3.613 


B 


600.00 


1.62 


18.03 ± 0.01 


243.9 ± 7.3 


Apr 2 02:31:56 


3.621 


V 


600.00 


1.61 


17.61 ±0.01 


322.6 ± 13.3 


Apr 2 02:44:17 


3.630 


I 


600.00 


1.61 


16.82 ±0.02 


449.2 ± 25.0 


A ~« o no.in.oo 
Apr z 02:19:38 


3.613 


T 
.1 


ODU.19 


1 CO 

l.OZ 


i a ni _i_ n i o 
16. Ul ± U.13 


613.6 ± ( (.1 


Apr 2 02:31:53 


3.621 


H 


560.18 


1.61 


15.40 ±0.09 


689.5 ± 68.2 








Epoch 3, St = 


34.9 min 






Apr 3 03:21:50 


4.656 


B 


600.00 


1.63 


18.62 ± 0.01 


142.1 ± 4.5 


Apr 3 03:34:01 


4.664 


V 


600.00 


1.66 


18.10 ±0.02 


204.1 ±9.0 


Apr 3 03:46:14 


4.673 


I 


600.00 


1.69 


17.34 ±0.03 


278.7 ± 16.4 


Apr 3 03:21:48 


4.656 


J 


480.15 


1.63 


16.56 ± 0.11 


368.7 ±41.3 


Apr 3 03:33:58 


4.664 


H 


560.20 


1.66 


15.88 ±0.03 


443.6 ± 23.4 








Epoch 4, St = 


34.9 min 






Apr 4 02:29:21 


5.620 


B 


600.00 


1.61 


18.54 ±0.02 


153.8 ± 5.5 


Apr 4 02:41:30 


5.628 


V 


600.00 


1.61 


18.08 ±0.01 


207.7 ±8.7 


Apr 4 02:53:44 


5.636 


I 


600.00 


1.61 


17.31 ±0.03 


285.8 ± 16.3 


Apr 4 02:29:18 


5.619 


J 


560.18 


1.61 


16.72 ± 0.17 


318.8 ±52.3 


Apr 4 02:41:28 


5.628 


H 


560.18 


1.61 


15.91 ±0.07 


431.9 ±33.6 








Epoch 5, St = 


33.2 min 






Apr 5 02:55:21 


6.638 


V 


600.00 


1.61 


18.59 ±0.03 


130.6 ±6.1 


Apr 5 03:07:36 


6.646 


I 


600.00 


1.63 


17.90 ± 0.03 


165.8 ±9.9 


Apr 5 02:44:50 


6.630 


J 


480.17 


1.61 


17.16 ±0.17 


212.5 ±35.4 


Apr 5 02:55:19 


6.638 


H 


560.19 


1.61 


16.43 ± 0.08 


266.0 ± 23.3 



Table 1 — Continued 



Date a At h Filter Exp. Time Airmass Mag. c Flux d 

UT day sec sec 2 piJy 

(1) (2) (3) (4) (5) (6) (7) 



Epoch 6, 5t = 34.6 min 



Apr 6 02:54:10 


7.637 


B 


600.00 




1.62 


19.57 ±0.02 


59.2 ±2.1 


Apr 6 03:06:13 


7.645 


V 


600.00 




1.63 


19.03 ±0.03 


86.6 ±4.0 


Apr 6 03:18:20 


7.654 


I 


600.00 




1.65 


18.21 ±0.04 


125.3 ±7.9 


Apr Uz.D4.Uo 


( .DO / 


T 

J 


OOU.l ( 




l.OZ 


1 / .04 zt U.Z1 


i oa v _i_ 07 a 
lOD. ( ± Z / .4 


Apr 6 03:06:11 


7.645 


H 


560.21 




1.63 


17.01 ±0.08 


156.9 ±13.3 








Enoch 7 St — 


34.7 








Apr 7 03:30:34 


8.662 


V 


600.00 




1.69 


19.20 ±0.03 


74.2 ± 3.5 


Apr ( Uo.4z.4o 


0.0 < U 


T 
1 


OUU.UU 




1 TO 
1. ( Z 


lo. 4o ± U.UO 




Apr 7 03:18:28 


8.654 


J 


560.18 




1.66 


17.81 ±0.22 


116.9 ±23.9 








Enoch 8 fit — 


34.6 


min 






Apr 8 01:48:29 


9.591 


B 


600.00 




1.63 


19.78 ±0.03 


48.9 ±2.0 


Apr 8 02:12:35 


9.608 


V 


600.00 




1.61 


19.21 ±0.04 


73.7 ±3.8 


Apr 8 02:00:33 


9.600 


I 


600.00 




1.62 


18.62 ±0.05 


85.4 ±6.1 


Apr 8 01:48:26 


9.591 


J 


1800.18 




1.63 


17.86 ±0.22 


111.4 ± 23.1 








Epoch 9, St = 


34.5 


min 






Apr 9 02:21:09 


10.614 


B 


600.00 




1.61 


20.01 ±0.04 


39.5 ±2.0 


Apr 9 02:45:13 


10.631 


V 


600.00 




1.62 


19.33 ±0.04 


65.7 ±3.4 


Apr 9 02:33:09 


10.622 


I 


600.00 




1.61 


18.75 ±0.06 


76.2 ±5.7 


Apr 9 02:25:49 


10.617 


J 


1680.18 




1.61 


18.00 ±0.21 


98.0 ± 19.7 








Epoch 10, St = 


36.4 


min 






Apr 10 03:05:21 


11.645 


B 


600.00 




1.65 


20.22 ±0.08 


32.7 ±2.5 


Apr 10 03:31:19 


11.663 


V 


600.00 




1.72 


19.59 ±0.07 


51.8 ±4.0 


Apr 10 03:17:38 


11.653 


I 


600.00 




1.68 


18.96 ±0.08 


62.4 ±5.7 








Epoch 11, St = 


36.0 


min 






Apr 11 02:47:16 


12.632 


B 


600.00 




1.63 


20.32 ±0.11 


29.8 ±3.2 


Apr 11 03:12:50 


12.650 


V 


600.00 




1.68 


19.76 ±0.09 


44.5 ± 4.0 
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Table 1 — Continued 



Date a 


At b 


Filter 


Exp. Time Airmass 


Mag. c 


Flux d 


UT 


day 




sec 




sec z 






(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


Apr 11 03:00:10 


12.641 


I 


600.00 




1.65 


19.23 ±0.12 


48.8 ±5.9 


Apr 11 02:47:14 


12.632 


J 


1800.19 




1.63 


18.44 ±0.26 


65.5 ± 15.9 






Epoch 12, St = 


33.4 


min 






Apr 15 02:58:11 


16.640 


J 


1920.15 




1.68 


18.94 ±0.32 


41.1 ± 12.3 






Epoch 13, St = 


36.2 


min 






Apr 21 23:36:46 


23.500 


V 


450.00 




1.90 


20.48 ±0.11 


22.9 ±2.5 


Apr 21 23:46:05 


23.506 


V 


450.00 




1.85 


20.50 ±0.12 


22.4 ±2.6 


Apr 21 23:55:43 


23.513 


I 


450.00 




1.80 


19.80 ±0.19 


28.9 ±5.3 


Apr 22 00:05:07 


23.519 


I 


450.00 




1.76 


19.42 ±0.14 


40.8 ±5.5 



Note. — These photometric observations have been grouped into thirteen epochs. The total duration of the epoch 
— from the time of the beginning of the first exposure until the ending time of the last exposure — is shown as St. 

a UT start time of the exposure. All observations were conducted in the year 2003. 

b Time since the GRB trigger. 

c Observed magnitude in the Landolt±Persson filter system, uncorrected for Galactic and host reddening. The 
errors given are statistical only. The systematic uncertainty in the zeropoint calibrations, to be added in quadrature 
with the statistical errors, are 5=0.03 mag, V=0.04 mag, 7=0.06 mag, J=0.05 mag, and i/=0.05mag. 

d Equivalent flux measurement at the effective wavelength of the respective filter (A e ff (B) = 4513.5 A, (V) = 5556.3 
A, (J) = 7935.0 A, (J) = 12440 A, (H) = 16528 A). As described in the text, the values for A off and zcropoints 
were determined self-consistently for the observed input spectrum. No correction due to reddening from dust has 
been applied. The errors reflect both the statistical and systematic magnitude zeropoint uncertainty in the flux 
measurement. No uncertainty in the filter flux zeropoint has been included. 



